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The problem of what we know, think we know, and think about the QCD coupling as is discussed. 
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1 The Ground gives rise to Measurements 

QCD has a split personality, "Perturbative QCD" and 
"Non-perturbative aspects of QCD" being routinely sep- 
arated by the organisers of HEP conferences. And so they 
are in our minds^. The microscopic dynamics of quarks 
and gluons is the QCD battleground; understanding the 
spectrum and interactions of hadrons is its ultimate ob- 
jective. 

The objective of this talk is less ambitious. My aim is 
to make you aware (if not convinced) of a possibility of a 
root that starts off with the QCD Lagrangian, employs a 
good old Dyson-Feynman field-theoretical staff of quark 
and gluon Green functions and may eventually lead to 
an understanding of colour confinement. To embark on 
such a quest one should believe in legitimacy of using the 
language of quarks and gluons down to small momen- 
tum scales, which implies understanding and describing 
the physics of confinement in terms of the standard QFT 
machinery, that is, essentially, perturbativelyB 

Is this programme crazy enough to have a chance to 
be correct? It seems it is. 

1.1 pQCD: a sketchy health report 

We shall start from a brief biased display of recent theo- 
retical advances on the perturbative frontier. 

1. Mass effects in heavy quark production xross sec- 
tions are now available at the NLO levelB 

2. A bunch of new state-of-the-art (2-loop -|- all- 
log-resummed) pQCD predictions have been de- 
rived to treat je± cross sectionsp hadrpproduction 
of heavy quarks □ r&nd prompt photonsP secondary 
hccivy quark pairsP the C-parameterl3 and broaden- 
ingLl distributions in e+e^ annihilation. Techniques 
are being developecLipr addressing next-to-next-to- 
leading order issuesEJ 



^ A reader willing to refresh his/her awareness of deep puzzles of 
the game is kindly adipised to consult the Proceedings of last year's 
HEP EPS conferencetJ 



Serious technical progress has been achieved in 
describing the High Energy Regime of scatter- 
ing cross sections of two small QCD objects, the 
BFKL Hemn (regretfully known as the "Hard 
Pomeron")llil This object should be responsible for 
a steep energy growth of production cross sections 
of 

HERA forward jets with pf ~ Q, 
widely separated in rapidity Tevatron jets, 
77 — > J/tp -f J/ip, and alike. 

The origin of the large-. NLO correction to the 
BFKL evolution kernel 113 is under scrutiny: how 
much of ii_is due to "kinematical" effects in the 
evolutionP^ running couplindifl angular ordering 
effects in space-like evolutionlla A physically moti- 
vated resummation of subleading kinematical and 
coUinear effects seems to gpaily improve the con- 
vergence of the FT analysisll^llj 

4. A unification of deep inelastic and diffractive phe- 
nomena is underway. It employs the notion of 
non- forward ("off-diagonal", "off- forward" ) [dou- 
ble] parton distributions which, on one hand, are 
related with various hadron form factors E3 and, 
on the other hand, can be accessed in hard in- 
teractions- such as deeply virtual Compton scat- 
tering ESE3 or hard diffractive electroproduction of 
(vector) mesonsEHl Of special interest are parton- 
helicity-jaensitive ("magnetic") non- forward distri- 
butions EiJ that do not contribute to the usual in- 
clusive DIS cross sections (structure functions) but 
participate, e.g., in determining the contribution of 
the ouark orbital angular momentum to the proton 
spinli3 see also t3 and references therein. 

5. Gluon radiation induced by propagation of a colour 
charge through a QCD medium (QGP, nuclear 
matter) attracts increasing attentionCJ 

6. An ideology and machinery for probing non- 
perturbative (np) effects with perturbative (pt) 
tools is being developed. 
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Table 1: World summary of measurements of as- Abbreviations: DIS = deep inelastic scattering; GLS-SR = Gross-Llewellyn-Smitli sum 
rules; Bj-SR = Bjorken sum rules; (N)NLO = (next-to)next-to-leading order perturbation theory; LGT = lattice gauge theory; resum. = 
resummed next-to-leading order. 







Q 










Aa 


,(M|) 






Process 


[GeV] 


a 


s{Q) 


a 




exp. 


theor. 


Theory 




DIS [pol. strct. fctn.] 


0.7 - 8 









1 20 + 0.010 


+0.004 
— 0.005 


+0.009 
—0.006 


NLO 




DIS [Bj-SR] 


1.58 





oyr + 0.062 
•J' — 0.081 





1 r,1 + 0.005 

^^-•^ - 0.009 






NNLO 


• 


DIS [GLS-SR] 


1.73 





r,qr + 0.092 

_ 0.073 





114+0.010 


+0.005 
— 0.006 


+0.009 
—0.010 


NNLO 


• 


T-decays 


1.78 





339 ± 0.021 





121 ±0.003 


0.001 


0.003 


NNLO 




DIS [v; F2 and F3] 


5.0 





215 ± 0.016 





119 ±0.005 


0.002 


0.004 


NLO 




DIS [fi; F2] 


7.1 





180 ±0.014 





113 ±0.005 


0.003 


0.004 


NLO 




DIS [HERA; F2] 


2-10 









120 ±0.010 


0.005 


0.009 


NLO 


• 


DIS [HERA; jets] 


10 - 100 









118 ±0.009 


0.003 


0.008 


NLO 




DIS [HERA; ev.shps.] 


7 - 100 









110 + 0.007 
— 0.006 


0.001 


+0.007 
—0.006 


NLO 




QQ states 


4.1 





223 ±0.009 





117 ±0.003 


0.000 


0.003 


LGT 


• 


T decays 


4.13 





220 ±0.027 





119 ±0.008 


0.001 


0.008 


NLO 


• 


e+e~ [tJhad] 


10.52 





20 ± 0.06 





1 on -1- 0.021 
J-OU fi n9Q 
— u.uzy 


+ 0.021 

n (190 


_ 


NNLO 




e+e~ [ev. shapes] 


22.0 





161 + 'fil 

— U.Uii 





1 -1- 0.009 
— u.uuo 


0.005 


+0.008 
— 0.003 


resum 




e''~e~ fcThadl 

^ ^ 1 ^-^ lltxU J 


34.0 





146 + 

— O.OZD 





1 90 -1- 0.021 
— 0.019 


+ 0.021 
— 0.019 




NLO 


• 


e~'~e~ [ev. shapes] 


35.0 





1 + 0.012 
^ — 0.007 





1 no -1- 0.008 
i-^iJ _ 0.006 


0.002 


+0.008 
—0.005 


resum 


• 


e~'~e~ [ev. shapes] 


44.0 





ion + 0.010 
^^'^ - 0.007 





1 90 -1- 0.008 
- 0.006 


0.003 


+0.007 
-0.005 


resum 




e~^e~ [ev. shapes] 


58.0 





132 ± 008 





123 ± 007 


0.003 


0.007 


resum 




pp ^ bbX 


20.0 





1 ir + 0.018 
— U.019 





113± 0.011 


+ 0.007 
— 0.006 


+ 0.008 
— 0.009 


NLO 




pp, pp 7X 


24.2 





1 07 + 0.017 

- 0.014 





111+ 0.012 

-^-^-^ - 0.008 


0.006 


+ 0.010 
- 0.005 


NLO 




(t(pp jets) 


30 - 500 









121 ±0.009 


0.001 


0.009 


NLO 


• 


e+c- [r(ZO ^ had.)] 


91.2 





122 ±0.005 





122 ±0.005 


0.004 


0.003 


NNLO 




c+c^ [ev. shapes] 


91.2 





122 ±0.006 





122 ±0.006 


0.001 


0.006 


resum 




e+e~ [ev. shapes] 


133.0 





111 ±0.008 





117 ±0.008 


0.004 


0.007 


resum 




e+e" [ev. shapes] 


161.0 





105 ±0.007 





114 ±0.008 


0.004 


0.007 


resum 




e+e~ [ev. shapes] 


172.0 





102 ±0.007 





111 ±0.008 


0.004 


0.007 


resum 


• 


e+e" [ev. shapes] 


183.0 





109 ±0.005 





121 ±0.006 


0.002 


0.006 


resum 


• 


e+e~ [ev. shapes] 


189.0 





109 ±0.006 





122 ±0.007 


0.003 


0.006 


resum 



The latter topic does not actually belong to the list be- 
cause its main objective lies beyond what is convention- 
ally considered to be the perturbative domain. 

1.2 Running QCD coupling - 1998 

A precise measurement of the coupling constant and ver- 
ification of asymptotic freedom remains a primary tar- 
get for experimental QCD studies. Not that chasing the 
third as{M^) digit would answer many (if any) a serious 
problepL Still Table 1 which summarises Us measure- 
ments Ej is very important for the theory. We should 
treat it with utmost respect as the major test of con- 
sistency of the tools that we employ to address various 
aspects of interactions involving hadrons. The results 
which appeared or were updated since summer 1997123 
(marked with a "•" ) include Us from 

1. the GLS sum rule, based on new GGFR v-N scat- 
tering datapa 



2. detailed^ and precise high-statistic ALEPH Ell and 
OPALlS studies of vector and axial- vector channels 
of hadronic r-decays, 

3. HI differential (2+1) jet rates at HERApi 

4. T-decays^ 



5. a CLEO determination from Er 



hadronic cross section measurement El 



52 GeV total 



6. reanalysis of JADE i?cm=35 and 44 GeV e+e" an- 
nihilation data to include ap_additional jet-shape 
variable — the C-parameterp3 

7. the most recent LEP update pf the ratio of hadronic 
to leptonic decay widthspj 

8. and finally, event shapes measured at-the highest 
LEP energies, E^rn=li>?> and 189 GevN 
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Figure 1: Compilation of as(Af|) measurements. 



Figure 2: Running QCD coupling-1998. 



The results on as{M^) fromE3Ho are still preliminary. 
Table 1 and Fig. 1 display as{M^) results evolved us^ 
ing the 4-loop Ms /^-function with 3-loop matching tB 
at quark pole masses Mf, =4.7 GeV and Mc =1.5 GeV. 
Fig. 2 collects values at the proper scales of individual 
experiments and spectacularly demonstrates asymptotic 
freedom. 

It is important to stress that "all meaningful sub- 
samples of results provide similar average values, and 
there is no significant systematic shift between any of 
those subsamples". Remarkedly, the "only e+e"" and 
"only DIS" subsamples peacefully coexist now, yielding 
a^Ml) =0.1210±0.0049 and as(M|) =0.1175±0.0061 
correspondingly. And so do "only Q <10 GeV" 
(0.1179±0.0043) and "only Q >30 GeV (0.1208±0.0058) 
measurements. 

The world average-s/alue of the ms coupling asiM^) 
is finally quoted to beEJ 

a,(M|) = 0.119 ±0.004. 

This result is based on 18 (of 27) individual measurement 
from the Table 1 with the errors Aa^ (M|) < 0.008. The 
overall uncertainty is derived by the "optimised correla- 
tion" methodE3 

1.3 as fresh from Lattice 

A lattice analysis of the-strong coupling was recently car- 
ried out by the grouped working on a QUADRICS QHl 



at Orsay. Calculations were performed on 16^ and 24^ 
lattices with a bare lattice coupling constant /3 = 6 (cor- 
responding to an inverse lattice spacing ~ 1.9 GeV). 
Calculations repeated for /3 = 6.2 (a~^ ~ 2.7 GeV) on 
a 24** lattice (which roughly embodies the same physi- 
cal volume as 16'* for /? = 6.0) yield consistent results. 
The primary objective of the study was to measure the 
coupling (Landau gauge, MOM subtraction scheme) "as 
defined in the text-books" . Namely, 

1. to measure 3- and 2-point gluon correlations, the 
Green functions 03(^1,^2,-^1 — P2) and G2{p) — 

2. truncate G3 to extract the vertex function 

3. and then construct the coupling at a symmetric 
Euclidean point >0, 

a,(fc2) = r3(pi,P2,P3)|p2^p|=p2^_fc2 • Z3/2(-fc2)_ 

Fig. 3 shows the result. The adjacent Fig. 4 displays the 
behaviour of a similar asymmetric correlator in which 
one of the three gluon momenta is set to zero. 

The fact that the lattice coupling shows the tendency 
to decrease at small momenta should not surprise us. 
Regardless of which particular kind of confinement we 
observe on the lattice, the correlators of coloured fields 
had better vanish at the origin. The gluon Green function 
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Figure 3: Running coupling on the lattice^ 



Figure 4: "Asymmetric" lattice coupling 







cannot have a pole at = 0, nor any similar singularity 
strong enough to propagate colour flux at large distances, 
j{R) ■ = O (1). Therefore, having chosen to truncate 
G3 with the text-book perturbative gluon propagators, 
we are bound to see Z{p^) trying to kill the massless 
gluon pole. 

It does not matter for our discussion how realistic 
these plots are. In particular, quarks were not included 
in the game. Even if they had been, there is always a 
question of how to incorporate light fermions with Comp- 
ton wave lengths comparable to (exceeding) the volume 
of the lattice world. This is an especially troubling ques- 
tion if light quarks do indeed play a crucial role irL the 
confinement of the world we live and experiment inB 

2 Measurements give rise to Assessments 

2.1 The origin and status of asik"^) 

Therefore, those who are not thoroughly aware 
of the disadvantages in the use of arms cannot 
he thoroughhL aware of the advantages in the 
use of armst^ 

We are used to the notion of a running QCD coupling. 
What is its origin and status? Formally, as is a param- 
eter of the perturbative expansion. From this point of 
view, its choice is, to a large extent, a matter of free will: 
it depends on how smart we want to be in organising the 
PT-series. It starts to make more sense, and becomes nat- 



ural, within a programme of trading a formal expansion 
parameter for a smarter object, a momentum-dependent 
as(fc^), which embodies some specific all-order radiative 
correction effects and is supposed to truly represent the 
interaction strength at a given momentum scale. 

A bunch of questions immediately comes to mind: 

• what high-order effects to embed? 

• what does ^Hruly represenf^ mean? 

• what is the characteristic momentum scale that the 
running coupling depends on? 

Renormalisability of the theory is known to be responsi- 
ble for the running of the coupling. The Renormalisation 
Group tells us: "scale the whole World, and the answer 
remains the same provided a has been changed accord- 
ingly" . Stretching the World is not at our disposal, how- 
ever. We may scale up/down external momenta but not, 
say, hadron masses. Therefore the running as rightfully 
appears in the realm of Hard Interactions characterised 
by a large momentum transfer, 3> m^, and in the 
description of observables insensitive to finite InfraRed 
parameters (particle masses m). 

The RG dictates how the renormalised coupling 
(constant) as(/i^) changes with the renormalisation scale 
/Ltij. It is casually said that in small-distance amplitudes 
characterised by a single large Euclidean virtual momen- 
tum scale Q 3> m, the formal dependence on the renor- 
malisation point, as(/i|j.), can be traded for as{Q^) thus 
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giving rise to the coupling (function) running with the 
physical momentum. 

A couple of cautious comments are due before we 
go further. Firstly, trading as{nj^) for Q!s(Q^) im- 
plies that hr/Q is the only momentum ratio that en- 
ters. This is true for renormalised off-shell ampli- 
tudes. These are however unphysical and consequently 
gauge-dependent objects; hence, the problem of peel- 
ing off gauge-dependent effects from the running cou- 
pling. Dealing with "physical" gauge-independent on- 
shell amplitudes has a different problem. They are in- 
frared sensitive and so have nasty as In Q /m (collinear) 
and a^lnQ/O (soft) high-order corrections. Therefore, 
one needs to treat collinear- and infrared-safe cross sec- 
tions instead of n-point amplitudes. 

Secondly, the way the running coupling emerges in 
Minkowskian observables is rather tricky. "On-shell" glu- 
ons are produced even from small distances 1/fc <C 1 fm 
with a coupling Q!<,(0^), in a manner of speaking, and 
not with as(fc^) (in a direct analogy with real photons 
from Z-peak being radiated with Q!em(0) ~ 1/137 rather 
than aem{Mz) — 1/128). What is instead determined by 
as(fc^) is an intensity of inclusive emission of a gluon sub- 
jet with the total transverse momentum k with respect 
to the radiating quark. 

At the two-loop level the as(/i|.)-dependence proves 
to be universal, i.e. independent of the choice of the pt- 
expansion parameter. This universality is inherited from 
the basic property of ultraviolet renormalisability of the 
theory. This part of the momentum variation of is 
determined by the first two terms in the /3-function, 

)-'.«„,,!),, 

namely, 

/3(a) = f3o + f3i^ + ... 

The further expansion terms /3„ with n > 2 remain 
scheme- (and even 5aM(7e)-dependent; in other words, ar- 
bitrary. Therefore, the large- momentum behaviour of the 
running coupling as{Q^) cannot be uniquely fixed be- 
yond two loops. The reason for that is pretty simple: 
only the first two loops are truly dominated by the uv 
region, that is by small-distance physics. 

There aren't many uv-dominated QCD diagrams. 
Logarithmically divergent radiative corrections to gluon 
and quark wave functions (self energies) and to gluon- 
gluon (quark-gluon) interaction vertices are involved in 
the renormalisation of the coupling constant a^- 

Consider for example a quark loop in the gluon self- 
energy. The one-loop radiative corrections contain the 
standard integral 

J ^ oc In Auv = oo =^ /3o ■ 



Hiding infinity under the carpet produces Po, the first 
coefficient in the pt /3-function expansion. In the next 
step we supply our quark loop with an additional inter- 
nal gluon. Now we have two independent loop-momenta 
to integrate over, qi and 92- Integration regions gi(2) ^ 
92(1) ^ Auv could have produced (InAuv)^ contribu- 
tions. These get suppressed by renormalising the inter- 
nal propagators and vertices at the one-loop level, the 
result being a single-logarithmic integral determined by 
the region qi 52 < Auv, 

--^as{fi%) oc as{fi%)lnAuv = 00 /3i . (1) 

This is how the usual story goes, order by order in pertur- 
bation theory. We can do better, however, by taking into 
consideration that the coupling in (^ runs with the in- 
ternal momentum. This means reorganising the pt series 
so as to incorporate into the two-loop diagram the higher 
order effects that result in substituting cks(/i^) —> as{q'^). 
By doing so we obtain a contribution which is still uv- 
divergent, though modified by the logarithmic decrease 
of the coupling at large momenta, 

— ^q;s((7^) cx Inln Auv = 00 =^ Pi ■ 

Renormalising it out gives rise to Starting from the 
third loop the situation however changes drastically: the 
uv-region is no longer dominant, and we get 




Pn>2 depend on the 
infrared physics! 



The notorious "nobody is perfect" applies to the above 
consideration as well. Strictly speaking, it is not known 
how to systematically refurbish, in all orders, the pt- 
series in terms of the running a^. So, this argument 
can be looked upon as yet annther example_Df a state- 
ment "correct but unproven"pl see howevercJ Still, the 
message it sends is clear: starting from the (next-to- 
next-leading) level, a purely perturbative treatment may 
become intrinsically ambiguous because of an intercon- 
nection between small and large distances. In particular, 
there may be no way of unambiguously defining the QCD 
coupling as, beyond the two loops, without solving the 
Theory in the infrared. 

2.2 Defining a^ beyond two loops 

Bearing in mind this warning one can still attempt to de- 
fine some PT expansion parameter as beyond two loops. 
To this end we may try different options. 

1. As long as higher terms in the /3-function are 
arbitrary anyway, why not simply set /3„>2 = 0. This is 
the so-called 't Hooft scheme which is perfectly fine but 
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for one thing: the Landau singularity — the fake infrared 
problem — obviously limits the use of as so defined, to 
sufRciently large momentum scales. 

2. One can design some sort of "optimisation prin- 
ciple" in order to fix ^2- The leading idea here is to 
minimise, one way or another, the effect of a (typically 
unknown, with few exceptions) three-loop pt correction 
to a given observable (Minimal Sensitivity, Fastest Ap- 
parent Convergence). Guided by sheer pragmatism (not 
necessarily a curse word) this approach often results in 
fixing P2 such as to force as(fc^) to develop a "spurious 
infrared fixed point" . 

A (curse) word spurious appears here because the 
3-loop /3-function and the corresponding coupling that 
emerge from such an analysis are observable-dependent. 
Linked to this is a well-grounded criticism based on non- 
transitivity between the couplings defined on the base of 
different observables. Moreover, it is legitimate to ask 
ourselves, whether it makes sense to optimise pt series 
which are known to be senseless, factorially divergent, 
being contaminated by infrared renormalons. 

Given all these reservations, it is still interesting to 
notice that such a "spuriously finite" coupling comes out 
close to the characteristic magnitude of the interaction 
strength in the infrared domain, the "couplant" <^^)) ~ 
0.2, which we shall discuss below. 

3. We may try to link as with some physical quantity 
in a "most, natural" way. To this family belong "effective 
charges" specific, physically motivated, perturbatiiie 
definitions aLa^ such as, e.g., the BLM^ CMW (MC)P 
or UraltsevEJ schemes fitted to serve gluon exchange po- 
tential, relativistic gluon radiation and non-relativistic 
heavy quark physics, respectively. The words "most nat- 
ural" don't deceive us. Whether we like it or not, there 
is no direct link between as and observables, the latter 
being hadronic observables, in a confining theory. 

4. ft is perfectly allowable to play with the /3- 
function (in particular, by adding to it non-analytic 
pieces) so as to make as(fc^) freeze (or even vanish) in the 
origin. Everything is allowed as long as we do not spoil 
the large-momentum behaviour of the running ^^(fc^) at 
the 1/lnfc and (f/lnfc)^ level. 

Rightfully rejecting the infrared "Landau singular- 
ity" in the coupling as being unphysical, or blaming in- 
frared renormalons for spoiling pt expansions, or dis- 
cussing possible links between as and hadronic observ- 
ables we stumble upon the very same problem: what is a 
true measure of interaction between inexistent objects? 
By hook or by crook, the problem is that of defining the 
theory in the infrared: the confinement problem. 

2.3 Does as at low momentum scales make sense? 
Possible responses are 



• Conservativei^No way. It cannot he, because this 
can't he, nevercB As long as there are no quarks and 
gluons in the physical spectrum of the theory, we 
cannot talk about the QCD interaction strength at 
distances above which colour-bearers get (mysteri- 
ously) confined. As a conservative position it is as 
impeccable as it is infertile. 

• Liberal: Let's try. Defining an infrared-finite cou- 
pling tts may serve as a tool fpt probing universal 
features of colour confinement^ {He who takes no 
risks drinks no champaign.) 

• Revolutionary: We must. Dyson-Feynman's is 
the only reliable Field-Theoretical language in our 
disposal. We have no other way but to employ 
quark and gluon Green functions down to A: = 
and search for an unusual, confinement, solution. 

Practically, God has generously supplied us with light 
quarks capable of delivering early screening, preventing 
colour JBelds from going haywire. As shown by V.N. 
Gribovp to achieve super- critical binding of light quarks, 
and thus to ensure the screening of any colour fields, it 
suffices to have the average magnitude of the coupling in 
the InfraRed region as large as 



> 



"crit 



0.f37. 



IR 



The interesting thing about this number (apart from be- 
ing easy to memorise) is that it is rather small. 

Pragmatically, if essential non-Abelian fields indeed 
grew really big, we would not even know how to properly 
define gluon degrees of freedom because of the notorious 
problem of "Gribov copies" e3 In Hamiltonian language, 
the Coulomb QCD interaction in the presence of trans- 
verse vacuum fields is described by the operator 



D-^^d'v + gs[Aj^,d^v], 



(2) 



reminiscent of the propagator of the Fadeev-Popov ghost 
in covariant gauges. In the second order in gs, the sta- 
tistical equal-time average of the product of two vacuum 
fields, (Aj^A^) cx ANc, takes over from the normal (uni- 
tary respecting) screening effect due to the splitting of 
the Coulomb gluon into "physical fields", namely two 
transverse gluons (—Nc/i) or a qq pair (— 2n^/3). 

This p hysi cal explanation of the anti-screening 
phenomenonpo had a dramatic continuation. The nor- 
mal magnitude of field fluctuations of spatial size L is 
LA-^{L) ^ 1. In the background of very large vacuum 
fields, LA^{L) ^ g^^ (so large that the non-linearity 
becomes essential at the classical level), the Coulomb 
(ghost) propagator (^) becomes singular. Technically, 
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this shows that we did not manage to properly define 
Lagrangian physical degrees of freedom, to divide out 
the volume of non-Abelian group transformations, to fix 
the gauge. 

A chilling perspective of Fadeev-Popov ghosts rising 
from the dead makes one wish to get away with a numer- 
ically small coupling (relatively small fields) as a unique 
(and maybe only) chance to keep things under control. 

Phenomenologically, there is no sign of strong colour 
fields: confinement appears to be soft and friendly. 
There are three aspects to this friendliness: 

1. pQCD works OK from very large scales down to 
Q ^ 2 GeV, for the phenomena where it should 
(hadronic r-decays being an extreme example); 

2. pQCD works down to (and below) Q ~ 1-2 GeV 
where it did not have to (as it does, in particular, 
for DIS structure functions at HERA); 

3. moreover, sometimes pQCD surprisingly works 
down to ... fc = NiLL. This happens, notably, in 
describing inclusive energy and angular distribu- 
tions of hadrons produced in hard processes, the 
phenomeDcm known as LPHD (local parton-hadron 
duality)ba'E2l 

What are the main lessons we have learned from experi- 
ment? 

1. We have got used to pQCD covering orders of mag- 
nitude in the basic hard cross sections. More than 
that, our toddler- wisdom of how to estimate np ef- 
fects is being certified. 

2. HERA tells us that proton is truly fragile. The 
transition from the jp physics to deep inelastic 7*p 
phenomena happens early, and it is sharp. The pro- 
ton isn't actually bound, if you take my meaning: 
a little scratch — and it is blown to pieces. 

3. Finally, when viewed globally, confinement is about 
renaming a flying-away quark into a flying-away 
pion rather than about forces pulling quarks back 
together. From the study of numerous string/drag 
effects in particle flows, in e"'"e~ and at the Teva- 
tron, we have learnt that even junky 200-300 MeV 
pions obediently follow the pattern of underlying 
colour fields. Whatever the ultimate solution of the 
confinement problem may be, it had better be gen- 
tle in transforming the quark-gluon Pointing- vector 
into the Pointing-vector of the final state hadrons. 



2.4 Some like it perturbative 

So, how to measure what we can't define? The idea is 
to look for deviations of inclusive quantities that char- 
acterise hard processes, from their respective perturba- 
tive predictions. For a ft- calculable observable such a 
deviation, due to genuine np (confinement) physics, is 
expected to be inverse proportional to a certain power of 
the hardness scale, 

Sa^ log'^g 

From within the PT-approach these (observable- 
dependent) powers can be inferred. Already at this stage 
we get a lot of valuable information. For example, 

• we would have no chance to see one of the most pre- 
cise determinations of as(M|) from the r-decay if 
confinement effects (which apriori could have been 
huge at as small a £jCale,as nir < 2 GeV) were not 
strongly suppressedEdS as (AQ(^j-,/m^)'^ <^ 1; 

• an understanding of the leading .nr power correc- 
tion to the Drell-Yan 7^-factorP§E3 l/Q^, calls 
for proper attention being paid to formulating the 
PT predictions for production of Drell-Yan lepton 
pairs, heavy flavours and high-invariant-mass jet 
pairs in hadronic collisions Liin order to avoid an 
artificial 1/Q contaminationfJ 

• discriminating jet-shape variables according to the 
envisaged power of the leading np contribution al- 
lows the selection of observables less affected by 
hadronisation and thus better suited for precise 
QCD tests such as measuring ag. 

As we shall discuss below, jet shapes typically contain 
l/Q NP contributions. Two examples of "cleaner" jet 
observables: the mean value of the three-jet resolution 
variable j/3, defined according to the Durham algorithm. 
which is hadronisatipn-scnsitive at the log Q/Q^ levelpZI 
or central moments^ ((^~(^))")j of nrost jet-shape 
observables 

However the story does not end here. Not only can 
one trigger the exponent of the confinement contribution 
to a given observable but, making further assumptions, 
one may quantify the absolute magnitude of the np effects 
by relating the deviations found for different observables. 

The programme can be looked upon as pushing for- 
ward the Sterman- Weinberg wisdom. PT-calculable ob- 
servables are Collinear-and-InfraRed-Safe (CIS) observ- 
ables, those which can be calculated in terms of quarks 
and gluons without encountering either coUinear (zero- 
mass quark, gluon) or soft (gluon) divergences. The pro- 
cedure is straightforward. 

''but not the broadenings, see later 
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1. Choose a CIS observable and get hold of the corre- 
sponding state-of-the-art PT-prediction;^ 

2. enjoy the beauty of the latter; 

3. observe that it has no sense (with as(fc^) running 
through the "Landau pole" ) ; 

4. force it to have one (kindly ask as(fc^) to behave); 

5. as far as the pt prediction produces now a definite 
answer, given an ir- finite as, see what you have 
done: quantify the ignorance about as in the low 
momentum range; 

6. do it again for other observables; 

7. verify that your ignorance is universal i.e. 
observable-independent ; 

8. eat your hat if it isn't and switch to another busi- 
ness. 

In recent years this programme has been carried out for 
a large set of practically interesting (majalities inthj, 
ing non-singlet DIS structure functimispZl'LBE^ e+i 
and DIS fragmentation functionsfj various jet shape 
characteristics (means and distributions in thrust, jet 
massea_£I-pa|Lanij3ter, jet hmadening) in e''"e~ annihi- 
lation SO'BBH and DISpO to hmvy hadron spec- 
tra froixL heavy-quark-initiated jets Lil and angular jet 
profilesllj First encouraging steps have been made to- 
wards revealing the pt-np interplay in hard small- a; 
phenonae|na such as high energy 7*7* scatterins-cross 
sectionjlS and non-singlet DIS structure functions 

Among interesting confinement-sick quantities still 
queueing for treatment are back-to-back energy-energy 
correlation and transverse momentum spectra of Drell- 
Yan pairs (which should exhibit weird non-integer Q- 
exponents), out-of-event-plane characteristics of two- 
and three-jet events (such as Tminor or oblateness), ac- 
companying Et hadron flows in DIS and hadronic colli- 
sions, and many others. 



3 Assessments give rise to Calculations 

How to quantify np corrections to hard observables? We 
shall sketch solutions to five major problems that arise 
on the way, namely 

1. How to disentangle power corrections coming from 
uv and m regions? 



'^These days, two-loop plus all-log rcsummcd with special care be- 
ing taken of the coupling running with the internal gluon momen- 
tum scale. 



2. How to split the magnitude of the power term into 
an observable-dependent vt- calculable factor and a 
universal np parameter? 

3. Is the latter really universal! What is the accuracy 
of the universality statement? 

4. How to merge pt and np contributions to the full 
answer? 

These being understood, the time comes to worry about 
your hat... 

3.1 Problem # 1: IR vs uv 

The trick of introducing small gluon mass m into Feyn- 
man diagrams canli 
momentum scales p^'"^^"^' 



robe contributions of small 
3 Operationally, the proce- 
dure is quite simple. Consider the pt correction to a 
given observable cr, given by one-loop Feynman diagrams 
with one additional gluon, real or virtual. Real and vir- 
tual contributions, and a^, are m-sensitive. Intro- 
ducing a finite mass into the Feynman propagator of the 
gluon regularises coUinear and infrared divergences pro- 
ducing 

2 

a,, a„ = ©(In^e) +0(lne)+0(l) , 

However, having chosen a CIS observable we have se- 
cured the cancellation of divergent terms in the physical 
answer, so that 

a — Or CTv — (1) 1 

thanks to the Bloch-Nordsieck theorem. By examining, 
how fast does the m-dependent contribution to a vanish 
in the e — > limit, we can quantify sensitivity of a given 
observable to large distances. To this end we havaJxi 
for contributions non-analytic in e in the origin[l3Lj 
while analytic corrections proportional to integer powers 
of rvi} IQ"^ come from the uv momentum region, k ^ 
in the Feynman integrals involving the modified gluon 
propagator, 



G{e) = 



1 



1 



-fc2 - iO 



iQ 



(3) 



For example, the one-loop total e+e annihilation cross 
section into hadrons has the following structure, 



1 



27r 



7^(e) , 7^(0) 



with 7?.(e) the "characteristic function" for the (Te+e- de- 
pending on the gluon mass. Setting e = we recover the 
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well-know first order pt correction uxjuBom — ols/t^- 
The function TZ is known, and so is its small-e behaviour, 



7^(e) 



11 



2e3 



Ine + O (e-^lne) . 



This shows that all the O (e) terms present in ct^ and 
Gy have cancelled in the sum, as well as the In^ e- and 
Ine-enhanced O (e^) corrections: non-analyticity starts 
at the O (e^) level only. This property can he easily 
inferred from a mysterious reciprocity relation e3 whose 
origin remains to be understood. 



7^(e) 



' 2 



The message that suelua powerful extension of the Bloch- 
Nordsieck wisdom e30 sends us is that the first gen- 
uine large-distance contribution to fJe+e- appears as a 
Q^^ power correction, corresponding to the ((-0?/;)^) vac- 
uum condfiusate of dimension six, in the standard OPE 
languagelI3'E2l 

We shall return to the issue of non-analyticity in m? 
a little later when we discuss the np trigger. But first, a 
confession is due: the last thing one would like to do is 
to make the gluon field massive. Violating sacred non- 
Abelian gauge invariance eventually results in breaking 
the uv renormalisability of the theory. 

Analytic coupling and "dispersive mass". Fortu- 
nately the situation is not so scary. What we actually 
need is not a real Lagrangian mass. It suffices to take 
into consideration the fact that the gluon, as any other 
quantum field, lives part-time in virtual states consisting 
of various parton configurations {qq, gg, etc). These vir- 
tual transitions appear as a renormalisation of the gluon 
field, 

1 Z(fc2) 



G(fc^) 



(4) 



which can be characterised in terms of the distribution 
over "virtual mass" m? (gluon spectral function) via a 
dispersion relation. It is this mass that makes the gluon 
"heavy". 

The factor Z has a simple physical meaning: it can 
be identified with the running coupling. 



(5) 



with /ifl an arbitrary uv-renormalisation point. Such 
an identification is straightforward in the Abelian the- 
ory: in QED, due to the Ward identity, radiative cor- 
rections to the wave function of a charged particle (elec- 
tron) and to the vertex cancel out {Zi = Z2)', what re- 
mains is the photon propagator, Z = Z^. In QCD the 



gluon which mediates interaction between "charges" is 
a charged (coloured) particle itself. Therefore its wave 
function renormalisation, unlike that of the photon, is 
gauge-dependent and has no direct meaning. What then 
enters into the definition of the physical coupling is the 
gauge- invariant combination of the gluon propagator and 
specific non-Abelian corrections to the interaction ver- 
tices, which are independent of the nature (colour repre- 
sentation) of the external charges. 



z = r(^^) • Z-, ■ r(^^) 



(6) 



Modulo this subtlety, in QCD we can still say that ex- 
changing a gluon with 4-momentum k brings in the prop- 
agator 



G{k) = 



-fc2 - iQ 



(7) 



(where we have dropped an irrelevant overall constant). 
Now comes the crucial assumption: we want (Q) to make 
sense in the entire complex /c^-plane. We know suffi- 
ciently well how as{—k'^) behaves in the deep Euclidean 
region, at large negative fc^; we know next to nothing 
about small-Zc^ region. However, whatever the function 
Us is, it had better respect causality. Therefore we sup- 
press the formal pt tachion (Landau singularity) and 
choose the "physical cut" alone, < fc^ < cx), as a sup- 
port for the dispersive relation. 



dm? q^ 
[m-^ + q^Y 



acff(m^) 



(8) 



Here Ofoff is the dispersive companion of the standard 
coupling, 

— ^aeff(Ai^) = -^Disc{a^(-/^2)| (9) 



din /i' 



The dispersion relation (^ can be formally inverted as 
the operator relation 



Q!eff(A' 



2 _ Sin(7r-P) 2 



-^«.(.^), r^,^—,. (10) 



It shows that in the pt region a^s and Us are pretty close: 
acff(fc^) = cts{k'^){l + {c?s))- One can design model ex- 
pressions for ois and the corresponding acS- What mat- 
ters in our discussion is that the latter gives us a direct 
hold on the large-distance contribution to the observable 
under study. We are ready for a "heavy gluon" . Substi- 
tuting (H) into the gluon propagator (Q) we can write 



as{-k^) 



dm 2 



din) 



- fc2 - iO 



Fcynman diagrams contain an integration over the inter- 
nal gluon momentum d'^/c. We conclude that the answer 
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can be represented as an Q;eff(r7T,^)-weighted integral over 
the "mass" of the (In m^-derivative of) the pt an- 
swer calculated by usual Feynman diagram techniques 
but with an as if massive gluon. Such a calculation pro- 
duces an m^-dependent characteristic function Ty spe- 
cific to a given observable V . In the first order in the 
running coupling is a function of e = m? /Q"^ and is 
given by the one-loop diagrams with a gluon of mass 
equal to the dispersive variable, < m? < oo. We get 

V{Q^,x)= ^a.^{m^)Me,xy, J'^--^, (11) 
m a In 6 

with X a set of relevant dimensionless variables. At this 
point ( pi] ) is no different from the usual pt answer. The 
CIS nature of V guarantees convergence of the inte- 
gration: vanishes as a power of e (e"'^) in the e — > 
(e oo) limit. Therefore the distribution Tv{e) has 
a maximum at e = Cvix) = C(l), and the integral is 
dominated by the large- momentum region ^ Q^. Ap- 
proximating acff(m^) ~ as{Q^) we easily reproduce the 
one-loop PT answer, 

aviQ^x) ~ as{Q^)TviO,x). (12) 

Using the observable-dependent position of the maxi- 
mum of the rn^-distribution as the scale for the coupling 
in ([ij), as{Cv{x)Q'^), does a better job since it minimises 
higher order effects. In this respect the dispersive tech- 
nology's closely related to the idea of "commensurate 
scales" N 

Renormalons. We could have tried better still by tak- 
ing into account pt corrections coming from logarithmic 
running of the coupling around the maximum. However 
there is an imminent danger around the corner. The de- 
ficiency of the PT series becomes apparent from the fact 
that the high-jXirder expansion coefficients Rk exhibit fac- 
torial growthE^ This is associated with both the uv and 
the IR integration regions in e. 

At the one-loop level we can equate acs with as and 
substitute for acsirn^) the standard geometric series 




= a,(Q2). 



This would give 

V{Q\x)-as{Q^)Tv{0,x) ~ asf^f^) Rk , 

k=l V ^ ^ 

with the expansion coefficients given by 

de f 1 \ ■ 
Rk= - fln-j Tv{e). (13) 



The ultraviolet contribution to Rk is estimated by inte- 
grating ( p^ ) over e > 1. If vanishes as e^' at large 
e, one finds for large k 

i?r-/°"*(ln^)' - {-q)-'k\. (14) 

This corresponds to an ultraviolet renormalon. Such an 
alternating series can be evaluated by Borel summation. 
This is because in the uv integration region the replace- 
ment aoff(TO^) — > Q!s(m^) is a reliable approximation and 
the contribution of this region can in fact be evaluated 
explicitly (e.g. numerically) without any expansion. 

The infrared contribution to Rk is estimated by in- 
tegrating over the region e < 1. The small-e behaviour 
of the characteristic function can be cast as 

Tv{e) ^ eP fv{\ne), (15) 

with / a polynomial at most quadratic in Ine. Both the 
leading power p and / depend on the observable under 
consideration. Substituting into ( p^ one finds 

i?iR = ^^^lnij eP f{\ne) ^ p-^ kl (16) 

We again find a factorial behaviour known as an in- 
frared renormalon. In this case however the coefficients 
are non- alternating and therefore the series is not Borel- 
summable. Attempts to ascribe meaning to such a series 
(which cannot even pretend at the status of asymptotic) 
give rise to unphysical complex contributions at the level 
of terms. This is generally interpreted as an in- 

trinsic uncertainty in the summation of the perturbative 
series. This problem is of a physical nature and can- 
not be resolved by formal mathematical manipulations 
alone. It requires genuinely new physical input to obtain 
a sensible answer. 

3.2 Problem # 2: np trigger 

To extract the np large-distance correction to a given 
observable we split the coupling into two components, 

a,(fc2)=aPT(fc2)+aNP(fc2)^ 

and use the latter as a trigger. It should be made clear 
that such a splitting is symbolic: it represents the cou- 
pling not in terms of two functions but rather of two pro- 
cedures. This ambiguity will be dealt with later on. For 
the time being let us just sketch these two procedures. 

• Having met a^^{k'^) under the integral we are 
advised to calculate it perturbatively, that is in 
terms of (not too long) a PT-series at the point 
fc^ ^ that our integral "sits" around. At the 
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same time we are supposed not to worry about our 
PT-coupling being potentially sick in the ir, momen- 
tum region. 

• On the contrary, integrals with a^^(fc^) are be- 
lieved to be determined by that very same IR re- 
gion, below some finite few-GeV scale. 



Since the integral of the type 
dp 



„NP,'j„2 



. fc2p = (few lOOs MeV)^^ (17) 



converges, it provides us with a dimensionful-coiistant 
(analogous to the ITEP vacuum condensateO for the 
OPE-controUed Euclidean quantities). These are the np 
parameters that will determine the magnitude of Q~^P- 
suppressed contributions to CIS observables. 

As long as a^^{k'^) falls fast enough, its dispersive 
companion a^g{m-^), according to should have van- 
ishing moments 



lim fcVP(fc2) = 



dm 2 



m2fa^ffP(m^) = 0, 



2\ 



at least for the first few integer values of p. That is 
why substituting a^^ for the full aeff in (0) singles out 
specifically non-analytic terms in the characteristic func- 
tion J-'v- The leading non-analytic term in the e ^ be- 
haviour of J^v determines then the power p of the 1 / Q^^ 
NP contribution. 

The usual view that the np component of the cou- 
pling decays fast at large momenta is rooted in the stan- 
dard ITEP-OPE ideologyka According to the ITEP pic- 
ture, NP physics is related to large-wavelength smooth 
gluon and quark fields. The presence of such fields in 
the vacuum does not affect the propagation of pt quanta 
with large Euclidean momenta, and as{k'^) in the uv do- 
main in particular. If those fields were non-singular at 
short distances, the quark and gluon propagators (and 
thus the running coupling) would be subject to expo- 
nentially small corrections only. Small-size instantons 
are believed to be the first (singular) non-trivial vacuum 
fields that disturb the propagation of quarks and gluons 
at the level of power corrections k^'^P with large expo- 
nents p > Pmax ~ /3o ~ 9. Phenonjipnological victories of 
QCD sum rules^ are impressiveHH This does not nec- 
essarily mean however that the true separation between 
the uv and ir physics is indeed as sharp and deep as it is 
implied by the ITEP picture. A viable alternative could 
be that the gluon Green function at large momentum 
fc, for example, bears the memory of the IR domain al- 
ready at the level of the relative fc"** term (the tail of the 

They [QCD sum rules] do not explain how the infrared soup is 
cooked^but taking this fact for granted, they skillfully utilise the 
recipec3 



IR-singular gluon polarisation operator) □ Such a heretic 
proposal does not necessarily undermine the notion of 
the gluon condensate- ( as G^). On the contrary, it can 
make it "calculable" In what follows we shall concen- 
trate on the practically most important smallest powers 
p = ^,1 and leave to future disputes the intriguing ques- 
tion of the depth of the separation between the uv and 
IR, phenomena. 

The non-analyticity of JF, necessary to generate np 
power correction, is typically of two kinds. In the first 
case we have an integer power accompanied by loga- 
rithm(s) of e which induces non-analyticity. This is the 
case of DIS structure functions, the Drell-Yan K-iactoT, 
the width of hadronic r-lepton decay and the total e+e~ 
annihilation cross section, etc. For example, for the va- 
lence DIS structure function one has 



hm Tms 

m — >0 



™2 g2 

a{x} ■ — In — - 



m 



Secondly, we may have a half-integer p. This is the 
case for many so-called jet-shape observables that charac- 
terise, in a CIS manner, the structure of final states pro- 
duced in hard processes. Thrust, invariant jet masses, 
G-parameter, jet broadening, energy-energy correlation 
(as well as a list of others not yet being discussed in 
the literature) belong to the P = | class. These quanti- 
ties shoulcLembody 1/Q power effects due to confinement 
physicsEJCa 



hm J^v — CLv ■ \ 
m->o \; 



Thrust. Historically the first among the jet shapes ad- 
dressed in this context has been the thrust, which mea- 
sures the "pencilness" of the final state system, 

T = max ' . 

2^^\Pi\ 

The direction n that maximises T is called the thrust 
axis. Two back-to-back particles with 4- momenta p and 
p (or clusters of particles with parallel momenta) pro- 
duced in the centre of mass of e+e" annihilation with 
cms energy Q correspond to T = 1 . Thrust deviates from 
unity for two reasons. One is pt gluon bremsstrahlung. 
T is obviously a CIS observable since neither coUinear 
parton splitting nor infinitely soft gluon radiation affect 
its value. Therefore the PT-component of the mean thrust 
is 

\(PT) 



(18) 



Another reason is pure hadronisation physics. pQCD 
radiation being switched off, two outgoing quarks are 
believed to produce two narrow jets of hadrons which 
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are uniformly distributed in rapidity and have limited 
transverse momenta with respect to the jet axis (Field- 
Feynman hot-dog, or a string). Within a simplified 
"tube model" with an exponential inclusive distribution 
of hadrons 



dN 
drjdk± 

we have 



i 

\Pzi I = 2 jdij dk_i 



dN 
dr^dk^ 

dN 
drjdk^ 



k±_ cosh 7] — 2fi sinh ry,„ = Q ; 
k±_ sinhry — 2/i(cosh77m — 1) , 



where, for the sake of simplicity, we have treated hadrons 
as massless. Constructing the ratio we obtain 

„ cosh??,„ — 1 2n ,^ ( 

T = — = 1 - + 0\ — 

sinh Q \Q'^ 

so that the departure of thrust from unity occurs at the 
1/Q level, with the characteristic hadronisation trans- 
verse momentum as a relevant scale. 



(l-T) 



(NP) 



2(fcj 



Q 



(19) 



Introducing finite hadron mass(es) does not change the 
result. It is from the study of hadronisation models that 
the 1/Q effects first came into focusEil 

The PT approach normally would not provide us with 
such a dimensionful parameter: gluon transverse mo- 
menta are broadly (logarithmically) distributed which re- 
sults in the mean {k±) cx asQ, in accord with ([l8|). 

However now we have a PT-handle on the large- 
distance physics, thanks to the "gluon-mass" trigger. 
Contribution to thrust from a single gluon with momen- 
tum k reads, in terms of light-cone (Sudakov) variables, 
k = ap + f3p + k_L , 

6(1 -T) = min{a,/3} . 

Calculating the characteristic function for (1 — T) one 
obtains 

J^r^ — [ —^dkl6{aPQ^ - fci - m^) • min{a,/3} 
TT J a p 



dk^, 



2Cf 



ttQ Jo yjk\ + VP?- 



Hence, 



dT 2Cf m 



-m 



dm? 



TT Q 

which is precisely the %/ m? non-analyticity we were ex- 
pecting. It is easy to see that the leading non-analyticity 



is due to the radiation of soft gluons at large angles, 
ko ^ k± ^ m <^ Q. 

Since soft gluon radiation has in fact a classical na- 
ture and its pattern is simple, as it is universal, the rela- 
tive magnitudes of the np contributions to CIS jet shapes 
appear to be simple as well. For example, for the mean 
jet shapes one obtains 



/TTxNP _ fly CV 
^ ' ^ Q ' 2tt 



dnn? 



NP/ 



(20) 



where the coefficients ay are simple numbers having a 
clear geometric origin. 

3.3 Problem # 3: Universality 

The very concept of the iR-finite coupling would make 
no sense without universality, i.e. if we had to introduce 
for each observable a private phenomenological parame- 
ter to fix the magnitude of the confinement contribution. 
Therefore it is natural thatJJan imiversality concept has 
been intensively argued forE^O'EJ 

It was soon recognised, however, that the tech- 
nologies based on the "massive gluon" trigger, whether 
the renormalon-motivated approach tracing the series of 
fermion bubbles or the dispersive approach, are intrinsi- 
cally ambiguous since there is no unique prescription for 
including finite-m^ effects into the definition of the shape 
variable. For example, a perfectly legitimate definition 
of thrust for a 3-parton system that consists of massless 
qq and a massive gluon, suggested by Beneke and Braun, 
has produced a result differing by a factor ~ 1.8 from the 
"conventional" (po|) The "naive estimate" that emerges 
as a result of substituting a "massive gluon" for the real 
final-state system of massless partons is fine for triggering 
the power of the NP-contribution but fails to predict its 
magnitude, the latter remaining prescription- dependent. 

Moreover, as has been already mentioned above, the 
running as{k'^) in Minkowskian observables emerges only 
at the inclusive level, as a result of an integration over 
positive virtualities of the gluon decaying into final-state 
offspring partons. Nason and Seymour rightfully ques- 
tioned the application of the inclusive treatment of gluon 
decays. They pointed out that jet shapes are not truly 
inclusive observables, since kinematics of the offspring 
matters for the V value. Therefore the configuration of 
offspring partons in the gluon decay may affect the value 
of the power term at next-to- leading level in as , which a 
priori is no longer a small parameter since the character- 
istic momentum scale is low. 

Both these problems called for analysis of the np- 
effects at the two-loop level. Such analysis has been 
performed, and the result came out unexpectedly sim- 
ple. It was shown that there exists a definite prescrip- 
tion for defining the "naive estimate" of the magnitude 
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of the power contribution, such that the two-loop effects 
of non-inclusiveness of jet shapes reduce to a universal^ 
observable-independent, renormalisation of the "naive" 
answer by the factor E3 

M ~ 1.76 (1.67) for = 3 (0). 

This is true for the NP-effects in the thrust, invariant jet 
mass, C-parameter and broadening distributions. The 
same factor (known as the Milan factor) also applies 
to the energy-energy correlation measure away from the 
back-to-back reeicin as well as to the linear jet-shape ob- 
servables in DISll3 

The universality of the Milan factor has three in- 
gredients. Firstly, it relies on the universality of soft 
radiation, the latter being responsible for confinement 
effects. Secondly, it obviously incorporates the concept 
of the universal the QCD interaction strength, all the 
way down to small momentum scales (which has been 
processed through the dispersive machinery). Finally, it 
stems from a certain geometric universality of the observ- 
ables under consideration, which includes their linearity, 
coUinear finiteness (convergence of rapidity integrals) and 
Lorentz invariance (independence of the gluan decay ma- 
trix element on the parent gluon rapidity) 

Strictly speaking, the accuracy of the Milan factor 
isn't great: the next loop would bring in the correction 



M 



true 



M 



2-loop 



l + O 



t)) 



with as entering at small scale. Hence, a ~ 20% uncer- 
tainty in the value of the A4 factor cannot be excluded. 
At the same time, the above ingredients of its universal- 
ity seem to be general enough as to ensure its validity 
even beyond two loops. 

3.4 Problem # 4- Merging 

We agreed to treat a^^ as a procedure rather than a 
function, and thus representing the NP-component of the 
answer in the form (2^) is not the end of the story. It 
is unsatisfactory in two respects. Psychologically, it does 
not satisfy our curiosity about interaction strength. More 
importantly, it remains symbolic since its PT-counterpart 
is given by a renormalon-sick series. We need to marry 
the PT and np components into a reasonable answer. To 
this end we first trade cteff back for the standard coupling 
by invoking the identity following from (0) , 



Jo 'm- TT Jo 



NP 



Then, we truncate the integration at some arbitrary fi- 
nite value fij which is large enough to neglect the np 

"^One can avoid introducing a^g in the firsfe-iplace and instead di- 
rectly exploit analytic properties of as (fc^)Ej 



interaction, a 
substituting 



NP 



(k) 0, k > ^i, and get rid of by 



dka'^^'ik') 



dkas{k^ 



dfca^^(fc^) 



The first integral quantifies the average interaction 
strength in the IR, via 



aoipi) 



1 

l^i Jo 



dk as{k^). 



The second contribution should be calculated perturba- 
tively, as a series in as- The series for this subtraction 
term is factorially divergent in high orders, as is the basic 
series for (V)^'^ . Introduction of the ir matching scale ni 
makes the full combined answer renormalon-free though. 
Finally, 

{V)^{Vf^{as)+avr, (21) 

where the universal np parameter V ^ 1/Q equals, to 
the second order in 



ACfM Hi 



Q 



/3o^ (ln^ + l + ^ 
2n \ Hi Pq 



(22) 



Here K is a known number depending on the scheme 
chosen for the pt expansion parameter as — as{Q^). A 
residual ^/-dependence, at the level of O [alfii/Q)j is 
the price for having a renormalon-free answer. In prin- 
ciple it can be reduced by continuing the PT-series, both 
for (V)^'^ and the subtraction in (p2|). The coefficients 
ay in ( ^l|) for mean thrust, C-parameter and invariant 
total and heavy-jet masses in e^e^ annihilation are 



V = 


l-T 


C 


m| 


Mil 


ay = 


2 




2 


1 



(23) 



The same NP-parametcr (^T]) that governs the leading 
confinement effect in the means affects the jet-shape 
distributions as well. In the "soft" kinematical region| 
Aqcd/Q <C V ^ Vmax it resieaJs itself as a shift in the 
pure perturbative spectrumpHo 



—iv) 

dV^ ' 



dV 



(24) 



The phenomenology of l/Q effects in jet-shape means 
and distributions suggests 



ZTT 







fc2 



^-^■kaTik 



0.2-0.25 GeV. 



*at the expense of introducing an additional phenomenological 
function, the shape function, it is possible to extend the range 
down to y = 0, asirecently demonstrated by G. Korchemsky in the 
test thrust caseEj 
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Leading power corrections to other observables are deter- 
mined by higher moments of a^^ . In particular, studies 
of power corrections to DIS structure functions al- 
low the quantification of the second mommt of the np- 
coupling, corresponding to p = 1 in (p^pjQ 



A2 = 



Cf 
27r 



■ea^^'ik') ^ 0.2 GeV" 



fc2 



Dedicated discussions of different aspects of what we 
know, think we know and think about as can be found 
in the literature. Among the mostjecent are: OPE, du- 
ality and the "physical coupling" Ip the simplest model 
for the IR- finite asEll a more sophisticated model which 
fits two known moments, P — and p = 1_ and re- 
spects 0PEp3 discussion of the causality issuepi instruc- 
tive lessons about OPE, m- finite ppupling and the Lan- 
dau pole from the 0{N) (T-modell£3 

4 Calculations give rise to Comparisons 

The phenomenology of power-suppressed contributions 
to jet shapes had a vibrant but somewhat troubled child- 
hood. Only thrust and C-parameter remained unaffected 
by theoretical misconceptions (some of which we are go- 
ing to sort out below). 

A ball-park value of ao — 0.5 was repeatedly emerg- 
ing from the analyses of jet shapes in p+e~ and DIS 
current jets in the Breit frameO'oEjO'EHl A typical re- 
sume of 1997 would run like "The concept of a 'univer- 
sal' Power Correction parameter ao in DIS ep scatter- 
ing and e'*"e~ annihilation is supported^ u3 However the 
same HI group was the first to complain about the ex- 
pectedin Q enhancement of the 1/Q correction to broad- 
ening which was inconsistent with experiment. The 
puzzle was recently put under scrutiny and clarified by 
the resurrected JADE collaboration. A detailed analy- 
sis presented at the Montpellier QCD conference by Pe- 
dro Movilla Fernandez showed that hadronisation effects 
in broadening not only shift the distribution to larger 
B values (as it is the case for 1 — T and C) but also 
squeeze it. As a consequence, fitting the data with a 
log Q-enhanced i3-independent shift yields inconsistent 
results for the total and wide-jet broadening distribu- 
tions (marked "old Bt" and "old Bw" in Fig. 5). They 
are also far away, in the parameter space, from the con- 
sistent pair of the thrust and C-Jiarameter distributions 
(solid T and C ellipses in Fig. 5)l£j Alarming experimen- 
tal conclusions "inconsistent results for jet broadening 
variables" and "effective coupling moments incompatible 
with each other" (read: there is no universality) made 
theoreticians jump on the train and revisit the problem 
from their side. 

It was soon recognised that one essential phe- 
nomenon was overlooked in the original NP-treatment of 



0.7 



0.6 



old B,. 



o 



i 0.5 



0.4 



\\Mh 



Bt 



0.3 



0.1 0.11 0.12 0.13 
as(Mz) 



Figure 5: 95% CL contours for jet shape means (dashed) and some 
distributions (solid). 



broadeninglxj namely an interplay between np- and pt- 
phenomena. To put it short, the effects produced by 
NP-radiation in the presence of normal PT-radiation are 
different from the effects of NP-radiation inferred from a 
pure first-order analysis, that is when the PT-radiation is 
"switched off" . 

NP-effects in the presence of PT-radiation. It is 

not for the first time that such "mistake" has been made. 
There is a whole list of confusions that came from the 
first-order analysis disregarding normal ever-present pt- 
gluons. 

The story of the heavy-jet mass is the simplest ex- 
ample. As you may remember from the previous dis- 
cussion, to trigger the NP-contribution we. are advised to 
add to the parton system a soft gluei^ a gluon with 

^ m ^ Aqcd- Let us do so at the Born level, that 
is add a gluer to the qq system as the third and only 
secondary parton. Constructing A'P of the quark-gluer 
system we will find a l/Q confinement contribution to 
the squared mass of the heavy jet, the one our giuer be- 
longs to. Meanwhile, the opposite lighter jet containing a 
lonely quark gets none of it. As a result the NP-correction 
to came out equal to that for thrust. 



(wrong) (25) 



In reality there are always normal pt gluons in the 
game which are responsible for the bulk of the jet mass: 
Mfj/Q'^ ~ a* > Ml/Q"^ > 6M^p. In these cir- 

cumstances it is not gluer's business to decide which of 
the jets is going to be heavier. Confinement effects are 
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instead shared equaUy, see (12 



2a 



M? 



(right) (26) 



It is worthwhile noticing that experimental analyses car- 
ried out before 1998 were based on the wrong expectation 
(p5|), which is obviously not the experimenters' fault 

Another still popular "mistake" is to expect that the 
NP-contributions to jet shapes can be suppressed by mea- 
suring higher moments, for example, ((1 — T)"). The one- 
gluon analysis indeed would formally produce for such an 
observable 

An 

((1 — T) )j^p Qj^^y ~ 

However what we are dealing with in reality instead is 

(i-r)" = [(i-r)pT + (i-T)Np]". 



Symbolically, 

(1-T)PT + 



Aqcd Y 

Q J 



A, 



QCD 



Q 



A 



QCD 

Q 



The leading 1/Q contribution is still here, reduced by 
the as{Q^) factor but far more important than the 1/Q" 
term. 

Another "mistake" of this sort brings us closer to the 
-B-issue. Consider the transverse momentum broadening 
of the current-fragmentation jet in DIS, that is the sum of 
moduli of transverse momenta of particles in the current 
jet. Adding a gluer to the Born (parton model) quark 
scattering picture we get three equal contributions to B: 
two contributions from the quark p which recoils against 
the gluer k emitted either in the initial (IS) or in the final 
(FS) state, \pj_\ = \k±\, and one contribution from the 
gluer itself when it belongs to FS. Taking into account 
the PT-radiation, however, the FS quark has already got 
a non-zero transverse momentum, p^"^ ~ Qs • Q, a sub- 
stantial amount compared to k± ^ Aqcd- In this envi- 
ronment the direct gluer's contribution is the only one 
to survive: the NP-recoil upon the quark gets degraded 
down to the l/Q^ effect after the azimuthal average is 
performed, 



{\p±\) = (wr - k^\) = p1 



„PT 



o 



Aqcd 



The true magnitude of the 1/Q contribution turns out to 
be a factor three smaller than that extracted from the 
one-gluon analysis. 

Nqw we are ready to address the squeezed broadening 



issuet 



STo the best of my knowledge the latest JADE analysis is the first 
one that properly included the Milan factor and fixed the Mh 
confusion e3 



The feature that 1 — T and C have in common is 
that the dominant NP-contribution to these and similar 
shapes is determined by radiation of gluers at large an- 
gles. This radiation is insensitive to the tiny mismatch, 
Qq = 0(as\ between the quark and thrust axis direc- 
tions which is due to pt gluon radiation. Therefore the 
quark momentum direction can be identified with the 
thrust axis. 

The broadening, on the contrary, accumulates con- 
tributions which do not depend on rapidity, so that the 
mismatch between the quark and the thrust axis starts 
to matter both in the B-means and distributions. 

If one naively assumes that the quark direction co- 
incides with that of the thrust axis, then B accumulates 
NP-contributions from gluers i with rapidities up to the 
kinematically allowed value ?/; < ?7inax — ln((5/A:ti). In 
this case one finds the shift in the i?-spcctrum to be loga- 
rithmically enhanced. 



A_B = OubV ■ In 



Qb ' 



(27) 



where ob ^ l(^) for the total (single-jet; wide-jet) 
broadening^ What was overlooked here is the fact that 
the uniform distribution in rji (defined with respect to 
the thrust axis) holds only for rapidities not exceed- 
ing |ln0q|. High-energy gluers with koi > ku/Qq are 
collinear to the quark direction rather than to that of 
the thrust axis and therefore do not contribute essen- 
tially to B. As a result, the NP-contribution to B comes 
out proportional to the quark rapidity. 



6Bl 



(NP) 



aiV ■ In 



(28) 



How does this affect NP-contributions to (B) and to the 
-B-distributions? The power correction to the mean sin- 
gle jet broadening (i?)^ is obtained by evaluating the 
perturbative average of SBi in 



{SB)\ 



(NP) 



aiV ■ ( In 



(29) 



At the Born level, the PT-distribution in the quark angle 
&q is singular at Qq = Q. In high orders this singularity is 
damped by the double-logarithmic Sudakov form factor. 
As a result, the NP-component of (B)-^ gets enhanced by 



ln± 



2y/CFas{Q) 



(30) 



For the mean wide jet broadening {B)y^ the result has 
the same structure with the replacement Cp ^ 2Cf due 
to the fact that now it is radiation off two jets which 
determines the 8g distribution. 
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The shift in the single jet (wide jet) broadening can 
be expressed as 

A,{B) ~ a,r-(\n^'^ , (31) 

where the average is performed over the perturbative dis- 
tribution in the quark angle Oq while keeping the value of 
B fixed. Since Qq is kinematically proportional to B, the 
log-enhancement of the shift in the i?-spcctrum becomes 

A^iB) 0^ aiV-ln^ (32) 

(with Bq a calculable function slowly dependent on 
a^lni?. Thus, the shift in the Bi (Bw) distribution 
becomes logarithmically dependent on B. 

The shift in the total two-jet broadening distribution 
At{B) has a somewhat more complicated B-depcndence. 
In the kinematical region where the multiplicity of gluon 
radiation is small, as In^ B <C 1, one of the two jets is re- 
sponsible for the whole PT-component of the event broad- 
ening, while the second is "empty". That "empty" jet 
contributes the most to the shift: in the absence of per- 
turbative radiation the direction of the quark momentum 
in this jet stays closer to the thrust axis. This results in 

At{B) ^ Ai(i?) + (i?)r' ^ (^Ini + i^^^^^ ■ 

In these circumstances the _B-dependence of the to- 
tal shift practically coincides with that of a single jet. 
In the opposite regime of well developed PT-radiation, 
as In^ B ^ I, the jets are forced to share B equally, and 
we have instead 

At{B) ~ 2 • Ai(S/2) ~ 2 • aiP In — . 

B 

It is the log B-enhancement of the np shift that makes 
the final distribution narrower than its PT-counterpart. 

On July 27^^ when this talk was being delivered, the 
origin of the "mistake" was already clear but the true 
answer still unknown. Therefore it wasn't risk-free to 
bet that the i?-distributions would eventually agree with 
respectable T and C. Fig. 5, which is preliminary, shows 
95% CL contours for as(M|) and ao as extracted from 
Bt distributions in the energy range 35-183 GeV, and 
also from the mean values of C, T, Mfj, Bt and B\y- 
Some care is to be taken in interpreting these results, 
since no accountJias been taken of the correlation of 
systematic errorsEj 

The universality of confinement effects has moved on 
to firmer ground, and with it, the concept of an ir- finite 
coupling. 



5 Comparisons give rise to Victories 

In our pursuit of confinement effects we were following 
the logic of the great Chinese warrior-philosopher Master 
Sun Tzu who said: 

The rules of the military are five: measurement, 
assessment, calculation, comparison and victory. 
The ground gives rise to measurements, measure- 
ments give rise to assessments, assessments give 
rise to calculations, calculations give rise ta^om- 
parisons, comparisons give rise to victoriesi3 
In our context, it would be a bit premature to talk about 
victories yet. It will suffice to have noticed that 

• QCD is alive and remains a rich source of technical 
and conceptual theoretical problems, of experimen- 
tal challenges and phenomenological amusement; 

• things are orderly in the "Hard Domain" , and 
as{Q'^) among them; 

• a crazy idea of probing perturhatively the "Soft Do- 
main" has gained ground; 

• the pQCD-motivated technology for triggering and 
quantifying, in a universal way, genuine confine- 
ment effects in hard observables is under construc- 
tion, 

• and that the effective large-distance inter- 
action strength inferred from these studies, 

0.14-0.17, turns out to be suffi- 
ciently small as to apply PT-language, at least 
semi-quantitatively, down to small momentum 
scales. 

According to my assessment, even if you 
have many more troops than others, how 
can that help you to victory?^ 
Heat is building up, and QCD is about to undergo a faith 
transition: we are getting ready to convince ourselves 
to talk about ^^quarks and gluons^^ down to, and into, 
the InfraRed. This is the core of tha^Gribov programme 
of attacking the colour confinement o It is interesting to 
remark that the average IR coupling, though numerically 
rather small, appears at the same time to be just large 
enough to activate the super-critical light-quark confine- 
ment mechanism. Is it a mere coincidence? Could be. 
To answer the question, the last Gribov works should be 
understood and developed (the write-up of the second 
paper concluding his two-decade QCD study remained 
unfinished). This won't be easy, given the complexity of 
the problem and, sadly, stone-solidness of our prejudice]^ 

''It was solid enough already 20 years ago when the semieal pa- 
per on what is known now as Gribov copies, G. horizonEI (459 
SPIRES citations, Dec. 98) was initially rejected by Nucl.Phys.B. 
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When on surrounded ground, plot. 
When on deadly ground, figh& 

It is worth to keep in mind that on the QCD ground, as 
we know it today, to fight isn't enough. To plot is often 
necessary. There is hardly a theorem around for which 
a proof could be concluded by the respectable QED = 
Quod Erat Demonstrandum (what had to be shown). 
The abbreviation QCD suits more. Quod Convenit 
Demonstrandum (what was agreed to be shown). A call 
for imagination and intuition, and guts to defend them, is 
what singles out the QCD island of the SM archipelago, 
the island where you never get bored. 
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